nature VJ>= — 

COMMUNICATIONS 



ARTICLE 

Received 29 Nov 2013 | Accepted 24 Jun 2014 | Published 1 Sep 2014 ^^^^^^^^^^3 OPEN 

A Cenozoic-style scenario for the end-Ordovician 
glaciation 

Jean-Frangois Ghienne\ Andre Desrochers^, Thijs R.A. Vandenbroucke-^, Aicha Achab^, Esther Asselin^, 
Marie-Pierre Dabard^, Claude Farley^, Alfredo Loi"^, Florentin Paris^, Steven Wickson^ & Jan Veizer^ 



The end-Ordovician was an enigmatic interval in the Phanerozoic, known for massive 
glaciation potentially at elevated CO2 levels, biogeochemical cycle disruptions recorded as 
large isotope anomalies and a devastating extinction event. Ice-sheet volumes claimed to be 
twice those of the Last Glacial Maximum paradoxically coincided with oceans as warm as 
today. Here we argue that some of these remarkable claims arise from undersampling of 
incomplete geological sections that led to apparent temporal correlations within the relatively 
coarse resolution capability of Palaeozoic biochronostratigraphy. We examine exceptionally 
complete sedimentary records from two, low and high, palaeolatitude settings. Their 
correlation framework reveals a Cenozoic-style scenario including three main glacial cycles 
and higher-order phenomena. This necessitates revision of mechanisms for the end- 
Ordovician events, as the first extinction is tied to an early phase of melting, not to initial 
cooling, and the largest 5^^C excursion occurs during final deglaciation, not at the glacial apex. 
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Shelf sedimentary architecture is controlled essentially by 
relative rates of base-level change and sediment supply^' . 
The base level reflects the interplay between tectonics 
(subsidence, volume change at mid-oceanic ridges) and the 
orbitally tuned, glacio-eustatically driven, sea-level change. The 
rate of the latter, at tens of metres per 10"*- 10^ years, is one to 
three orders of magnitude greater than the tectonically driven 
sea-level change, at tens of metres per million years or less. The 
critical issue for analysing the stratigraphic record, therefore, is 
the correct assignment of depositional units to their appropriate 
temporal hierarchy alongside a given sea-level curve. Another 
consideration is the temporal significance of the observed or 
suspected hiatuses. Any stratigraphic record of ancient shelf 
deposits, and their isotopic or palaeontological proxies, inevitably 
samples only the discontinuous segments of a given sea-level 
curve^, which often are below the relatively coarse resolution 
correlation potential of Palaeozoic biostratigraphy"*. Regardless, 
shelf deposits are the principal record that we have for pre- 
Mesozoic glaciations and they must therefore serve as 
stratigraphic archives for glacially driven events, providing: 
subsidence was active; water depths at the onset of glaciation 
were moderately deep; and sediment supply was adjusted to 
subsidence rates. These preconditions are essential for the 
maintenance of significant water depths during glaciation, as 
any rapid shallowing would pre-empt the registration of 
subsequent glacio-eustatic events. 

The end-Ordovician witnessed one of the three largest 
Phanerozoic glaciations with the development of continental- 
scale ice sheets^"^. This climatic event was postulated to have 
been initiated by massive weathering of fresh volcanic rocks^, 
tectonics and related plate motions^" high cosmic ray flux 
impacting cloud albedo or by a combination of the above^^. The 
glaciation apparently coincided with highly^"* or moderately^ ^'^^ 
elevated CO2 levels, with large isotopic excursions (C, S, O, N, 
Nd), and with a major double-phased biological extinction 
Interpretations based on far-field, low palaeolatitude sequences, 
resulted in 'coup-de-theatre' scenarios that have tied the two 
phases of extinction to the onset and termination of a single 
glaciation ^^'^^'^^ Yet the high palaeolatitude near-field sequences 
contain up to five glacial cycles that can be tentatively correlated 
across the Gondwanan glaciated platforms^'^^. The low 
palaeolatitude archives must therefore represent a more 
complex scenario ^^'^^'^^"^^ than that of a single, major glacial 
event. If so, a multiorder climate signal with a hierarchy of cycles, 
a Cenozoic-type 'business as usual' scenario, is a more likely 
alternative than a large singular event. Such linkage of eustatic, 
biological and isotopic records to the climatically forced 
development of an ice sheet can only be contemplated within a 
framework of high- resolution sequence stratigraphy that 
integrates alio-, chemo- or biostratigraphic markers. 

Here, we present such a framework, based on the recognition 
of genetic stratigraphic sequences (GSSs) and intervening erosion 
surfaces (see Methods). This framework, driven by glacio-eustatic 
cycles tied to the evolution of polar continental- scale ice sheets 
over west Gondwana^, enables the correlation of eustatic cycles at 
a level that is beyond the resolution capability of most absolute 
dating methods and of biozones, the latter typically of Myr 
duration^. A Cenozoic-style scenario including three main glacial 
cycles and higher- order phenomena necessitates the revision of 
the end-Ordovician, glaciation -related sequence of events. 

Results 

Palaeolatitude sequence stratigraphic frameworks. We intro- 
duce sequence stratigraphic correlation frameworks for two 
superbly exposed and exceptionally well-developed latest 



Ordovician successions (Fig. 1), the Anti- Atlas of 
Morocco^'^^'^^ and Anticosti Island in Canada^^'^^. Both offer 
sections, on a 100 -km scale, from the basin edge to the axis of 
active sedimentary depocentres (Fig. 2 and Supplementary Figs 1 
and 2). Relative to the end-Ordovician ice-sheet centre (present- 
day north-central Africa), they provide a near-field (Anti-Atlas, 
siliciclastic platform) and a far-field (Anticosti Island, mixed 
carbonate and siliciclastic) stratigraphic record. These two 
successions, up to 300 and 100 m thick, respectively, were 
deposited in basins with notable subsidence rates and significant 
(ca. 100 m) initial Katian water depths, enabling the development 
of comprehensive archives of the latest Ordovician glaciation 
(Supplementary Fig. 3). On the basis of average shelf- depositional 
rates within the overall Late Ordovician context^'^ , and on 
comparison with analogous late Cenozoic shelf stratigraphies^'^^, 
such thick successions are considered to be long-term 
»100skyr archives. In both areas, the end-Ordovician 
comprises three genetic low- order stratigraphic sequences 
(GSSs) of highest significance that, in turn, encompass a 
number of higher- order GSSs of intermediate and low 
significance (Fig. 2; see the Methods). 

The intercontinental correlation of these two successions is 
made possible due to the recognition of marker intervals. Earlier 
palaeontological studies already bracketed the stratigraphic 
interval that contains the well-known end-Ordovician extinction 
events (Supplementary Figs 1 and 2). In both sections, the first 
extinction event is situated around the conventional Katian- 
Hirnantian boundary, which in our record is penecontempora- 
neous with the major bounding surface that separates the two 
lower low-order GSSs. The related 'maximum flooding interval' — 
rather than the maximum flooding surfaces (MFS) that cannot be 
strictly synchronous at the global scale — is our first marker. 
It correlates with the brief 'pre-Hirnantian deepening' event 
identified in western Laurentia^^. The second marker based on 
allostratigraphy relies on the signature of the end-Ordovician 
glacial climax. In the Anti-Atlas, it demonstrably correlates with 
the stratigraphic interval bounded by glacial erosion surfaces and 
includes widespread glacial (subglacial, glaciomarine . . . ) deposits 
in the basinal succession (the glacial interval in Fig. 2). Note that 
coeval strata are often absent in basin edge successions (Hajguig 
Wadi log in Fig. 2). In the Anticosti Island succession, the 
signature of the glacial climax (lowest sea levels) is ascribed to the 
prominent erosional unconformity at the base of the Laframboise 
Mb. The result is a severe erosional truncation of the studied 
interval (Fig. 2 and Supplementary Fig. 2). Regional correlation of 
low- and high-order GSSs between these two markers is indeed 
intriguing (Fig. 2). Moreover, other subordinate erosional 
unconformities at the basin edge of the Anti-Atlas succession 
have their counterparts in the Anticosti Island succession. For the 
highest- order GSSs, at least partially related to local processes, 
such correlations are less reliable. 



The end-Ordovician glacial tempos. Within the context of 
glaciation, where eustasy is expected to control shallow shelf 
sequences^, our findings strongly suggest that the two 
independent regional scale frameworks and their correlation are 
robust and that the correspondence of the low- and high-order 
GSS records from dissimilar tectonic and environmental settings 
arises from glacio-eustatically fluctuating sea levels, the latter a 
consequence of waxing and waning of the western Gondwana ice 
sheet. We interpret the three low-order GSSs to be the signature 
of the three extensive glacial cycles (Latest Ordovician Glacial 
Cycles, LOGCs 1-3; Figs 2 and 3 and Supplementary Note). Note 
that our provisional numbering refers to the latest Ordovician, 
understood to informally include the highest Katian and the 
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Figure 1 | Geological settings, (a) Position of the two study areas on a palaeogeographic reconstruction (modified from Achab and Paris^^). In the Late 
Ordovician, the Anti-Atlas succession was part of the North Gondwana shelf, an epicratonic domain at high palaeolatitudes (green star). It included an actively 
subsiding depositional trough, which was free of ice until the middle Hirnantian, but glaciated during the glacial maximum in the upper Hirnantian. In contrast, 
the Anticosti Island succession accumulated in a foreland basin that developed at low palaeolatitudes during closure of the lapetus Ocean at the onset of the 
laconic orogeny (orange star), (b) On Anticosti Island, cliffs and marine platforms provide a superbly exposed latest Ordovician to early Silurian mixed 
siliciclastic-carbonate shelf succession. At Pointe Laframboise, alternating limestone (inner ramp) and softer shale (outer ramp) intervals across the shoreline 
platform clearly reflect the sequence stratigraphic correlation scheme inferred from sedimentary facies analysis, (c) Extensive outcrops of the Anti-Atlas 
Palaeozoic succession offer a superb record of predominantly shallow-marine Cambrian to Carboniferous sequences^^'^^. Near Tazzarine, a thick (up to 350 m) 
succession of alternating sandstones (storm- to tide-dominated, subordinate glaciation-related deposits) and offshore shales allows a compelling sequence 
stratigraphic correlation framework to be established for the latest Katian and Hirnantian (see Fig. 2). Coloured triangles represent the three identified Late 
Ordovician Glacial Cycles (L0GC1-3) that comprise a lower, orange/yellow triangle (regressive system tracts, RSTs) and an upper, purple triangle (transgressive 
system tracts, TSTs). LOGO are bounded by dotted lines, underlining MPS. Other surfaces: dashed lines, sharp-based surfaces and their correlatives within RSTs; 
solid lines, maximum regressive surfaces; wavy lines, subaerial erosion (Anticosti) or basal glacial erosion surface (Anti-Atlas). 



Hirnantian. This is in agreement with the views that ice sheets 
were extant already before the latest Katian^'^^'^^'^^. Our first 
glacial cycle spans the upper Katian (LOGC 1), the second 
(LOGC 2) includes the uppermost Katian strata and most of the 
lower to middle part of Hirnantian and the third (LOGC 3) 
commences in the upper Hirnantian and ends in the lowermost 
Silurian. In our view, corroborated recently by Nd isotope 
studies^^, the end- Ordovician glaciation could not have been 
restricted to a single short-lived glacial event, as earlier believed. 

The minimum depositional time for the entire LOGC 1-3 
succession is in excess of the Hirnantian duration (~ 1.4 ±0.2 
Myr^^), the latter encompassing about 60-90% of the LOGC 2 
and some 40-80% of the LOGC 3. Assuming that all LOGCs 
are of about equal durations, a single LOGC corresponds to a 
0.7-1.6 Myr time span. The embedded higher-frequency 
multiorder event stratigraphy is typical of orbitally controlled 
climatic oscillations that lead to recurring ice- sheet growth 
stages^"*, in agreement with the modelling results of Hermann 
et al^^ Note however that in contrast to the well-known, strongly 
asymmetric and shorter-term. Pleistocene glacial cycles^^, our 
LOGCs show no abrupt deglaciations. They have a symmetric 
distribution of the high-order GSSs, as evident from the stacking 
patterns within the low-order regressive to transgressive system 
tracts (TSTs). Long-lasting interglacials are expressed as 
condensed, maximum flooding stratigraphic intervals^ that 
account for significant portions of the overall duration of the 
studied time span (Fig. 3). Despite of some similarities to 



Quaternary glaciations^'^'^^'^^'^^, the durations and internal 
organization of LOGCs argue for dissimilar glacial tempos and 
forcings. These Ordovician features and tempos more closely 
resemble the Oligocene climate patterns that were driven by a 
high -amplitude obliquity modulation at 1.2 Myr frequency^"*, 
resulting in a limited number of short-lived ice- sheet growth 
phases, our high-order GSSs, centred around the obliquity 
nodes^^'^^. Such high-frequency signals may hold some 
similarities to the metre- scale cycles described from other low- 
latitude areas and attributed to ^200 (ref. 40) or 40-130 kyr 
(ref. 32) frequency oscillations. 

Assuming the analogy with the Oligocene climate is valid, we 
hypothesize that an orbital forcing responding to the amplitude 
modulation of the obliquity typifies glacial climate systems at 
relatively high C02atm levels. In such a scenario, the ice-sheet 
inception, driven by ice-albedo feedbacks, may have resulted from 
a dearth of exceptionally warm rather than a ubiquity of 
exceptionally cool summers^^. 



Discussion 

Our sequence stratigraphic framework allows Hirnantian excur- 
sion(s) and extinction(s) to be revisited. The large positive carbon 
isotope excursions of the Palaeozoic, such as the Hirnantian 
Isotopic Carbon Excursion, HICE in LOGC 3 (refs 17,18,41-43), 
are often used as chronostratigraphic markers, albeit with no 
consensus model for their existence. Yet, the notion that the 
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Figure 2 | Sequence stratigraphic correlation frameworks in the Anti-Atlas and Anticosti Island successions. The two frameworks were established 
separately and later correlated (note the different vertical scales: 50 m and 100 unit bars for the Anti-Atlas and Anticosti Island successions, respectively). 
Both show three superposed low-order GSSs that connprise alternating RSTand TSTs symbolized by coloured, orange/yellow and purple triangles, 
respectively. These low-order GSSs include higher-order GSSs (high- and highest-order GSSs) of intermediate and low significance, symbolized by smaller 
triangles (white and grey triangles, respectively). Apparent intercontinental correlation of the two sequence stratigraphic frameworks for 
low- and high-order ranks confidently supports the proposition that the low-order GSSs are signatures of three successive LOGCs, which together 
encompass a 2-4 Myr latest Katian to lowermost Silurian time interval. Because the high- and highest-order GSSs relate to the relative rates of 
base-level changes (Supplementary Fig. 3), rather than to a priori assumptions about sequence duration, no temporal significance can be attributed at this 
stage. We retain here lithostratigraphic names because the related boundaries crosscut timelines, for example, the Katian/Hirnantian boundary (see 
Supplementary Figs 1 and 2). Lf: Laframboise Member. Black arrows point to stratigraphic interval with significant faunal turnovers. Encircled numbers 1-6 
refer to the numbering of glacial erosion surfaces (see Supplementary Fig. 1). 



S Qarb signal of shelf carbonates is a direct reflection of the 
S^^Cdic of the globally uniform open ocean is clearly open to 
debate^^'^^ (Supplementary Discussion). Note that the magnitude 
and occurrences of such enrichments depends on localized 
settings (for example, epeiric versus open ocean aquafacies^^'^^) 
and is therefore related to depositional facies and not 
straightforwardly to a global signal. For example, the 5^^Ccarb 
on the modern Bahamas Bank is considerably more positive than 
that of the open ocean"^^. 

In addition, our revised chronology questions the paradigm 
of temporal relationships that link the position of the end- 
Ordovician glacial cycles, their tempos and biochemical 
events^^'^^"^^'^^. The first issue that arises is the identification 
and temporal range of the HICE itself. If it is understood as 
coeval with the large +4%o isotopic excursion, it has to be 
confined to a restricted time interval of a single high-order GSS 
within the end- Ordovician glaciation (Fig. 3), as posited by the 
Anticosti case study. If, on the other hand, understood as a ^^C 
signal that commences in the latest Katian and ends in the latest 
Hirnantian, our results (Fig. 3 and Supplementary Table 1) show 

4 



^^C enrichments in at least three stratigraphic positions, suggesting 
that HICE combines several excursions, thus challenging its validity 
as a high-resolution chronostratigraphic marker. 

The Anticosti S^^Qarb curve (Fig. 3) includes two main 
isotopic events. First, it is the well-known excursion in the 
Laframboise Mb. ( + 4%o) that is disconnected from a rising limb 
in the underlying strata by a major unconformity that we relate to 
the glacial maximum and to subaerial erosion in LOGC 3. 
Second, there is an earlier asymmetric excursion ( + 2%o) with its 
descending limb that is spanning the lower and middle parts of 
the Ellis Bay Fm. (LOG 2). There is also a lesser enrichment in the 
uppermost Vaureal Fm., associated with LOGC 1, which may 
form a third, subordinate excursion. Other putative (<l%o) 
excursions, while present, are minor and difficult to interpret. 
This multi-peak isotope pattern at Anticosti Island questions the 
views of strictly synchronous signals, despite observations that a 
number of Hirnantian records worldwide — and potentially 
similar 'wiggles' in the carbon isotope record elsewhere — 
contain positive 5^^C spikes that appear isochronous ^^'^"^'^^ 
within the correlation capabilities of the Palaeozoic bio- and/or 
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Figure 3 | Perspective of sequence stratigraphy. Temporal correspondence between documented (chitinozoa^^'^^'^^'^^) or essentially inferred 
(graptolite) biostratigraphies and (a), the Anticosti Island succession with its related isotopic signal and faunal turnovers and (b), the Anti-Atlas succession 
with interpreted low-frequency sea-level changes and ice-sheet occurrences. A cycle hierarchy is developed that distinguishes L0GC1-3 
(low-order, high-significance Late Ordovician Glacial Cycles represented by both coloured triangles and the thick, pale blue curve) from high-order cycles 
(thin, dark blue curve). LOGCs are bounded by major MFS (dotted lines). Blue shading highlights time intervals specifically characterized, or thought to be 
characterized by ice-sheet development stages. The ice-sheet development increased from the late Katian to the late Hirnantian, as suggested by 
glacioeustatic trends. The dashed blue curve is a representation of the early Silurian eustatic background. Black, dashed lines are the inferred Katian to 
Hirnantian and Hirnantian to Llandovery boundaries, (c) Representation of a potential time calibration is based on astronomical forcings dominated by 
1.2 Myr amplitude modulation of obliquity cycles"^^ (see text for details). By analogy with the Cenozoic, the composite artificial curve was constructed by 
mixing high-frequency orbital cycles ('ETP' for eccentricity-tilt-climatic precession-^-^) and here it is shown only to illustrate the distortion in the 
stratigraphic record. It results in condensed transgressive and overdeveloped lowstand intervals, relative to a linear timescale. The high- and highest-order 
glacial cycles likely correspond to such orbitally forced, high-frequency climatic oscillations. In contrast, during the long interglacials orbital forcing did not 
result in ice-sheet development and they have therefore a poorly differentiated record. The end-Ordovician includes short glaciation intervals with 
cumulative duration of perhaps a few hundred thousand years. The embedded isotopic and biological signals show up to three discrete isotopic events and 
two faunal turnovers (oblique-line shading), from the highest Katian to uppermost Hirnantian. The Hirnantian isotopic carbon excursion (HICE) is not 
restricted to the excursion associated with LOGO 3 at the top of the Ellis Bay Fm. The dashed pink curve is a representation of the Katian istopoic 
background. 



chronostratigraphy^^. An apparent single peak may represent 
only disjointed parts (Fig. 4) of a hypothetically complete S^^Qarb 
curve for just one of several repetitious LOGCs, or a composite 
signal from an artificially stacked section. 

Whatever the temporal extent of HICE, our sequence 
stratigraphic framework warrants reconsideration of the pub- 
lished 'cause-and-effect' scenarios for its origin. The rising limb of 
the ^^C excursion at the base of the Ellis Bay Fm. (Fig. 3) is 
associated with a highstand that follows the LOGC 1-2 transition. 



while its descending limb spans several high-frequency glacio- 
eustatic cycles within the late regressive to TSTs of the LOGC 2. 
In this case, there is therefore no apparent connection between 
eustasy and the 6^^Ccarb curve. The simplest explanation is to see 
the LOGC 2 isotopic signal as that of regional epeiric water 
masses with their distinctive variations in 5^^C (ref. 44). In 
contrast, the subsequent, exceptionally high -amplitude excursion 
is within the TST of LOGC 3, and is associated with a drastic 
basin- scale change of facies caused by transition from glacial to 
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Figure 4 | A detailed interpretation of the far-field LOGC 3 stratigraphic and isotopic record, (a) The main lithostratigraphic units on Anticosti Island, 
including the Laframboise Member, are shown with their representative depositional facies and related 6^-^C curve. They are separated by shoreline 
ravinement unconformities^^. The sequence stratigraphic interpretation differentiates low-order/high-significance regressive (orange and yellow triangle) 
and transgressive (purple triangle) system tracts. The corresponding first-order unconformity coincides with the base of the Laframboise Member. 
High-order cycles represented by white triangles are present in the LOGC 3 TST, which commences with the base of the Laframboise Member, 
(b) The same sequence in the 'linear' timescale perspective of this succession of depositional events. It includes relatively long depositional hiatuses 
(oblique-line shading). At Anticosti, the high-frequency glacio-eustatic sea level changes, similar to those recorded in the near-field glacial record of 
Morocco, are represented by unconformities coeval with the glacial maxima. One recorded interglacial event (Laframboise Mb.) is interpreted here as a 
single high-order GSS within the larger TST of LOGC 3. (c) An alternative view of the isotopic excursion, which includes the stratigraphic hiatuses. The 6^-^C 
record captures only disjointed segments of the isotope signal. In particular, the 6^-^C curve does not include values from the time interval that corresponds 
to the Hirnantian glacial climax. We suspect that the trend from background levels to the maximum in fact combines an initial rise that predates the glacial 
climax, the associated hiatus and the subsequent maximum that postdates the glacial climax. This maximum is developed mainly within the reefal 
limestones constituting the highstand facies of a particular high-order GSS. 



warmer climates (reefs of the Laframboise Mb.). At a higher 
resolution, the excursion appears to be confined to the highstand 
of a high-order GSS (Fig. 4), thus peaking at times of rising sea 
levels associated with deglaciation. This coincidence is opposite to 
the postulated lowstand conditions that are essential in the 
weathering' scenario^ ^ and the model can be discounted as a 
potential explanation. The 'productivity"^^ and related 'circulation 
pattern' explanations"*^'"*^ could perhaps offer plausible alter- 
natives, providing it can be demonstrated that the isotopic 
excursion is not facies dependent. Our highstand nadir of isotope 
excursion can then be consistent with the scenario that invokes 
carbon storage in the deeper parts of the shelf"^^, albeit 
constrained— because of it high amplitude— to basinal, not 
global, scales (see box model in Supplementary Discussion and 
Supplementary Tables 2 and 3). In such a context, the particular 
highstand conditions favourable for the development of carbon 
excursions may arise at distinct locations during any high- order 



GSSs. If so, it is the short duration of contiguous high-order 
GSSs that give the impression of a synchronous, worldwide 
phenomenon during the LOGC 3 transgressive trends. For minor 
excursions, such as those in LOGC 1 or in the uppermost (below 
the unconformity) Ellis Bay Fm., we contend that our present-day 
knowledge of carbon isotope systematics does not permit unique 
diagnostics of causative factors and scenarios. We therefore desist 
from their interpretation. 

In summary, providing our sections represent sufficiently 
comprehensive archives of the latest Ordovician development, we 
dispute the apparent association of each LOGC with an individual 
isotopic excursion. At higher resolution, the relationship with sea- 
level history differs from case to case, indicating that it is not a 
unique forcing but likely a combination of processes that is 
involved in *^C enrichment^*. 

Similar reasoning may suggest that 'pulses' in patterns of the 
end-Ordovician biological extinction result from telescoping of 
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segments from the stratigraphic record, as must be the case for 
hiatus -dominated successions^^. The two phases of the Late 
Ordovician mass extinction that were documented worldwide in 
earUer studies ^^'^^ are however confirmed also by our results on 
the Anticosti Island (Fig. 3 and Supplementary Fig. 2). Our 
correlation framework moreover indicates that these turnovers 
are relatively long-lasting time intervals that encompass several 
glacio-eustatic fluctuations of the high or highest GSSs in 
the Anti- Atlas record. Whether these turnovers originate from 
evolution specifically related to LOGC developments, or whether 
they only mirror a succession of stacked, quasi- instantaneous, 
pulses is beyond resolution of our data set, regardless of their 
potential combination within protracted, global events^ The 
ensemble of our sections studied represents only a fraction of 
affected palaeohabitats and biota, and, as explained, sections that 
contain complementary palaeontological data cannot be readily 
correlated into our framework; our analyses thus likely 
undersamples the full biotic dynamics through this interval. 
Yet, our juxtaposition of extinction phases to glacial development 
suggests a more nuanced scenario than previously advocated 
(Fig. 3). 

The older turnover, which has classically been associated 
with the onset of the Hirnantian glaciation at the base of the 
N. extraordinarius Biozone, spans an interval that includes LOGC 
1 deglaciation and the early LOGC 2 highstand. This turnover 
corresponds therefore mostly to the first major interglacial period. 
The models that are based on processes linked to glacial onset, 
such as the shrinkage of biotic ecospace, temperature decline 
or development (or the loss) of anoxia during falling sea 
level^^'^^'^ '^^ are thus not compatible with this revised 
scenario. Instead, processes linked to deglaciation dynamics (for 
example, amplification of meltwater fluxes that enhance ocean 
stratification), or flooding of the shelves by relatively deep anoxic 
waters ^^'^^'^^ appear to be more likely scenarios for this first 
turnover, but they are not applicable, as previously envisioned, 
for the second turnover. The second turnover that we recognize 
in LOGC 3 is an event traditionally assigned to the lower part of 
the N. persculptus Biozone. This extinction/ recovery pattern 
affects mostly macrofauna in the Anticosti Island succession 
(Supplementary Fig. 2). The phytoplankton crisis, on the other 
hand, commenced beneath the regional unconformity (Fig. 4), 
that is before glacial climax of LOGC 3 (refs 42,52), suggesting 
that the second end-Ordovician faunal turnover may have been 
initiated already during the late Hirnantian ice-sheet waxing, thus 
casting doubt on a unique causative linkage that would have been 
confined to final deglaciation. Note nevertheless that the 
ubiquitous existence of worldwide hiatuses at that time makes 
any interpretation tentative. 

While we appreciate the merits of a sophisticated model-driven 
approach, and welcome the impetus derived from it, the insights 
arising here from the application of basic geological methods 
underline the need for detailed understanding of the rock record 
as well^^. In this contribution, temporal relationships of near-field 
and far-field records for the end-Ordovician glaciation are 
considered within a high-resolution, multiorder correlation 
framework that reflects a Cenozoic- style hierarchy of glacio- 
eustatically driven oscillations consisting of three main cycles and 
superimposed higher- order subcycles. This interpretation 
questions earlier views that were based on lower resolution data 
sets for a simple latest Katian decline in sea level followed by its 
rise in the upper Hirnantian. An oversimplified sedimentary 
succession likely incorporates significant hiatuses and represents 
only a partial record of the entire time interval^^'^^. Frequently, 
such stratigraphic sections correspond to vertically juxtaposed, 
unrelated parts of glacial cycles, resulting in biased timing of 
biochemical signals relative to the glacial tempos. On the basis of 



our framework, we anticipate that the most easily captured phases 
in the worldwide end-Ordovician development should reflect the 
initial waxing stage and potentially the immediately ensuing 
flooding event (LOGC 1), followed by later reflooding of the 
shelves at the end of LOGC 3. In more comprehensive 
successions, the maximum flooding interval at the LOGC 2-3 
transition (mid-Hirnantian transgression^) will likely yield a 
decipherable signature^^. 

The orbitally controlled depositional record of a glacial interval 
will mostly be underrepresented in proven or suspected hiatuses. 
These may originate from nondeposition, subaerial or subglacial 
unconformities, transgressive post-glacial ravinement processes, 
mass movements or from erosion by bottom currents, the latter 
being particularly effective for the deeper parts of shelf basins. 
Due to the lack of Palaeozoic deep-sea records, the absolute 
timing and calibration of the Ordovician glaciation may remain 
enigmatic^ ^. We envision therefore that future progress in 
understanding the temporal, spatial and causal evolution of the 
Late Ordovician environmental record will have to rely on high- 
resolution methods that capture multiorder sequence stratigraphy 
and related proxies along depositional profiles on a regional scale. 
Our advances using these methods include: rejection of the earlier 
cause- and- effect scenarios for HICE(s), as these no longer fit with 
the revised context of glacial/ glacio-eustatic development; the 
suggestion that low-order LOGCs likely represent the 1.2 Myr 
long obliquity cycles that modulated ice-sheet dynamics, similar 
to scenarios proposed for the Oligocene; and the insight that the 
first Hirnantian extinction pulse, contrary to earlier studies, was 
linked to an intervening melting phase, not to the initial cooling 
phase of the end-Ordovician. 

Methods 

GSSs. Sequence stratigraphic correlation frameworks are based on visual corre- 
lations of marker beds along continuous exposures at the 10-30 km scale (Fig. 2) 
and on refined, regional scale, chitinozoan-based biostratigraphies for northern 
Gondwana^^"^^ and eastern Laurentia^^'^^. This results in correlations that are 
noticeably different from lithostratigraphic schemes (Supplementary Figs 1 and 2). 
MFS and a variety of erosion surfaces have been delineated in the field. The MFS 
coincide with deeper, usually condensed, depositional conditions and serve as 
bounding surfaces for GSSs^^ ideally including a lower regressive system tract 
(RST) and an upper transgressive system tract (TST). Erosional surfaces 
correspond to glacial erosion surfaces (Morocco); subaerial unconformities 
reworked by transgressive ravinement processes (SR-U sensu Embry, 2009 (ref. 62); 
Anticosti); or sharp -based erosional surfaces punctuating regressive facies trends 
and ascribed in most cases to regressive surfaces of marine erosion (Anti- Atlas and 
Anticosti). We favour GSSs over Trangressive/Regressive^^ (T-R), or depositional 
sequences^ because their bounding surfaces (MFS) better approximate late 
deglaciation conditions and thus appropriately bracket glacial cycles. In this 
scheme, a post-glacial highstand of an interglacial is represented by deposits that 
constitute the lower part of the subjacent sequence. 

Sequence hierarchy. Stratigraphic surfaces have been assigned to a hierarchy of 
GSSs. The significance of facies shifts and/or their penetration into the basin are 
used as criteria to assess the relative magnitude of base-level falls in successive, 
multiorder sequences. It results in a data-driven hierarchy^^, different from a 
frequency- related scheme based on a priori assumptions about durations of 
sequences. The highest-order (low significance) GSSs display limited facies shifts, 
both in the basin and at basin edge. More significant are the high-order genetic 
sequences, which comprise several highest-order GSSs and/or include abrupt facies 
shifts associated with coeval, or at least suspected, erosion surfaces at basin edge. 
The low- order sequences (highest significance) are made up of a suite of high -order 
sequences, the stacking pattern of which defines long-term RST and TST. They are 
bounded by the major MFS associated with severe condensation (for example, 
phosphogenesis in the Anti- Atlas). They include in their most regressive part (late 
RST or early TST) one or several important erosional surfaces such as shoreline 
ravinement unconformities, or glacial erosion surfaces in the upper Hirnantian in 
Morocco, which expand toward basinal areas. This approach is often not 
appropriate for maximum flooding intervals characterized by relatively deep 
depositional conditions, where facies shifts are poorly deciphered. Here, an 
alternative, frequency- related, hierarchy is frequently applied^. 

Base-level falls associated with glacial erosional surfaces are recognized on the 
basis of: their basin ward extent at regional scale^^ (Supplementary Fig. 1); the 
development/absence of well- organized subglacial shear zones that indicate fully 



NATURE COMMUNICATIONS | 5:4485 | DOI: 10.1038/ncomms5485 | www.nature.com/naturecommunications 

© 2014 Macmillan Publishers Limited. All rights reserved. 



7 



ARTICLE 



NATURE COMMUNICATIONS | DPI: 10.1038/ncomms5485 



subglacial/marginal ice fronts^'^. Maximum erosional depths are not considered to 
be a measure of the significance of a glacial surface. We are aware that such 
estimates reflect glacial extents rather than true ice- sheet volumes, but they do have 
significance when dealing with high- and low-order GSSs. 
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